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Abstract —Linear detectors such as zero forcing (ZF) or 
minimum mean square error (MMSE) are imperative for 
large/massive MIMO systems for both the downlink and nplink 
scenarios. However these linear detectors require matrix inver¬ 
sion which is computationally expensive for such huge systems. 
In this paper, we assert that calculating an exact inverse is not 
necessary to find the ZF/MMSE sointion and an approximate 
inverse would yield a similar performance. This is possible if the 
quantized solution caicniated nsing the approximate inverse is 
same as the one calculated using the exact inverse. We quantify 
the amonnt of approximation that can be tolerated for this 
to happen. Motivated by this, we propose to use the existing 
iterative methods for obtaining low complexity approximate 
inverses. We show that, after a sufficient number of iterations, 
the inverse using iterative methods can provide a simiiar error 
performance. In addition, we also show that the advantage of 
using an approximate inverse is not limited to linear detectors 
but can be extended to non linear detectors such as sphere 
decoders (SD). An approximate inverse can be used for any 
SD that requires matrix inversion. We prove that application of 
approximate inverse leads to a smaller radius, which in turn 
reduces the search space leading to reduction in complexity. 
Numerical results corroborate our claim that using approximate 
matrix inversion reduces decoding complexity in iarge/massive 
MIMO systems with no loss in error performance. 

I. Introduction 

With growing demand for high throughput, Mutiple-Input- 
Multiple-Output (MIMO) systems with large/massive number 
of antennas are expected to become an indispensable part of 
fifth generation wireless technology Q, El. It employs a 
large number of antennas at the base station (of the order 
of hundreds) that operate to serve relatively fewer users. 
However, we know that as the number of antennas grow, the 
complexity of detection algorithms increases 0. Thus, there is 
need for techniques which, while exploiting the extra degrees 
of freedom, are able to decode the transmitted signal efficiently 
in terms of error performance and complexity. 

In the literature. Zero Forcing (ZF) and Minimum Mean 
Square Error (MMSE) have commonly been used as precoders 
in a massive MIMO downlink a, Q and as decoders in 
a massive MIMO uplink. Even the complex decoders for 
uplink transmission also require the computation of ZF/MMSE 
solution. For example, neighborhood search based algorithms 
a, El or sparsity based detectors 0, in use such linear 
detectors for initialization. Calculating a ZF or an MMSE 
solution requires inversion of a matrix. However, finding an 


inverse is computationally expensive, especially when large 
number of antennas are employed. 

In this paper, we argue that an approximate matrix inverse 
suffices for finding a ZF/MMSE solution. In other words, 
usage of an approximate inverse does not compromise the 
quality of a ZF/MMSE solution. Since the solution obtained 
using linear detectors anyway needs to be quantized, it is clear 
that there is a scope for using an approximate inverse as long as 
the quantized solution remains unchanged. We derive bounds 
on the approximation such that the ZF/MMSE solutions from 
the exact and approximate inverses are same. Further, we show 
that the advantages of using an approximate inversion are not 
limited to linear detectors. Thus, a class of Sphere Decoding 
(SD) algorithms ifTol require the ZF solution for computing the 
Bahai Radius (BR) lUIl, II2, consequently requiring matrix 
inversion. Hence, one can think of utilizing an approximate 
matrix inverse even in complex decoding schemes like SD. 

In this work, we propose the application of an approximate 
inverse to compute the BR for usage in SD. The approximate 
inverse has two advantages. Firstly, it reduces the complexity 
of matrix inversion. But secondly, and more importantly, 
we prove that it results in a smaller BR. This is a bigger 
advantage as complexity of decoding in such SD algorithms is 
largely governed by the choice of BR. Simulations results for 
large/massive MIMO systems corroborate that the proposed 
SD provides a low complexity solution with no loss in error 
performance. 

II. System Model 

Consider a massive MIMO downlink with N transmit 
antennas at the base station and K users, each with a single 
receive antenna. Such a system can be represented by 

yd = HdSd + nd, (1) 

where = Wx^, W is the linear precoder such as ZF or 
MMSE and is the N dimensional signal vector transmitted 
from the base station. Each element in x^ is drawn from a set 
H, all entries of which belong to an M-QAM constellation, 
with average symbol energy Eg. represents the K x N 
i.i.d. channel matrix with zero mean and unit variance and 
is an i.i.d. zero mean Gaussian noise vector with dimension 
K xl and variance Nq. The i-th entry of the vector y^, yi^d, 
is the signal intended for the i-th user, for / = 1, 2, • • ■ ,K. 


Similarly, in the case of uplink, the system can be repre¬ 
sented by 


inequalities are satisfied by the error matrix E (a sufficient 
condition) 


y„ = H„x„ -f n„, (2) 

where x„ is the K dimensional transmitted signal vector 
whose i-th entry is the symbol transmitted by the i-th user, for 
i = 1,2, ■ ■ ■ ,K. Again, each element in x„ is drawn from the 
set fl, with average symbol energy Eg. Similarly, H„ is the 
N X K i.i.d. channel matrix with each coefficient having zero 
mean and unit variance. The noise vector n„ is i.i.d. x 1 
Gaussian with each element having zero mean and variance 
Nq, and Yu is the N dimensional received signal vector at the 
base station. This results in KEs/Nq Signal-to-Noise Ratio 
(SNR) at each receive antenna. 


III. A Linear Detector using Approximate Matrix 
Inverse 

Linear detectors such as ZF and MMSE are useful for both 
the uplink and downlink (as a precoder) in massive MIMO 
systems. The expressions for these detectors can be expressed 
as 


XZF 

Xmmse 


(H^H+^lK)-'H^y 

-t^S 


( 3 ) 

( 4 ) 


where [-J quantization operator to the set il and H is the 
N X K channel matrix. Quantization allows us to use an 
approximate inverse instead of exact inverse while giving the 
same ZF/MMSE solution. Since the operations are similar in 
both the uplink and downlink scenarios, we consider only the 
uplink scenario for the analysis. Eor notational simplicity, we 
have removed the subscripts here onwards. 

Let us define the error in the approximation of the inverse of 
a matrix C as E = C — C“^, where C = H^H is a AT x AT 
matrix that needs to be inverted and C is its approximate 
inverse. Also, define g = H^y. 


A. A Bound on the Acceptable Error in the Matrix Inverse 

We will consider an approximate matrix inverse good if the 
ZE solution calculated through it is equal to that calculated 
through the exact inverse. In this section, we evaluate a bound 
on the error which can be tolerated in the computation of an 
approximate matrix inverse. 

Eor the ZE solutions calculated using the exact and approx¬ 
imate matrix inverses to be equal, the following equality must 
be satisfied 


argmin||x- C-^H^ylP 
=> argmin||x — C“^g|p 
argmin||x — C“^g|p 

xGH 


argmin||x — CH^yjp 

argmin||x — Cg|p 
x^f2 

argmin||x — C“^g — Eg|p. 

( 5 ) 


Let the solution of the L.H.S. of (|5]l be Xzf and let z = 
XzF — C“^g. Therefore, (|5]l will be satisfied if the following 


, K 

^min cvA/ \ ' 7-1 \ ^min 

2 ^ ~ < 2 
i=i 

I^min N dmin 

2 ^ . Ejjdj) ^ 2 ’ ^ 

i=i 

Vi = 1, 2, ■ • • , fc and j = 1,2, ■ ■ ■ , k, where Zi is the i- 
th element of z, E^j is the (i,j)-th element of matrix E, 
dmin is the smallest distance between any two points in the 
constellation, and 3? and 3 denote the real and imaginary parts 
respectively. After combining the K equations in (|6]l and (|7]l 
and taking expectations on all sides, we have 

^lK<E[K(z-Eg)]<^l^, 

^lrc<E[S(z-Eg)]<^l^. 

where Ik is a AT x 1 vector with all entries as ones. 

Here E(zi) = 0, for all i = 1,2, ■ ■ ■ , N, because for a 
given transmitted vector x, Xzf can take any point in the 
constellation around x due to randomly and independently 
distributed noise and hence the expectation of difference 
between the two quantities would be zero. Therefore 

^l;c<E[5R(Eg)]<^l;, 
which, after substituting for E and g, yields 

< E[5?((Cfc - C-i)(H"Hx + H^n))] < 

^ < E[sR(Sx)] < ^ 1 k , (8) 

where we define S = I —CC as the residual matrix. Similarly, 
^1 k<E[3(Sx)]<^1k (9) 

Hence, if ([8]l and (|9l) are together satisfied by S for a given 
transmitted vector x, the ZE solutions through approximate 
and exact inverses would be equal. Next, we discuss some 
low complexity approximate matrix inversion methods which 
can be used to find ZE solution accurately. 

B. Low Complexity Iterative Methods for Computing Approx¬ 
imate Matrix Inverses 

Several low complexity iterative methods for finding the 
inverse of a matrix have been proposed in ifTSll . lfT4l . Let Cj^ht 
the approximate inverse and be the residual matrix after k 
iterations. The order of the iterative method is p if the residuals 
after k and k 1 iterations satisfy = S^. Eor e.g., in 
a third order method, approximate matrix is calculated in the 
following manner ifTsll 

Cfe+i =Cfe(3I-CCfe(3I-CCfc)), 


( 10 ) 




















where I is the identity matrix. Here, we note that S^+i = S^. 
Similarly, a seventh order iterative method is defined as 

Cfc+i = Cfe(7I + CCfe(21I + CCfe(35I + CCfc(35I 
+ CCfc(21I + CCfe(7I + CCfe)))))), (11) 

and here, we have S^+i = S^. 

In our simulations, we use Newton’s iterative method for 
finding approximate matrix inverse which has low latency, 
low complexity Eol and is also easy to implement m. The 
approximate inverse is updated in each iteration according to 

Cfe+i = (2I-CfeC)Cfe. (12) 

Here, S^+i = S|, revealing quadratic convergence. Increasing 
the number of iterations increases accuracy, but also increases 
the number of operations required and hence affects complex¬ 
ity, resulting in a trade-off between performance and efficiency. 

Initial matrix Cq needs to be chosen with care as it decides 
the number of iterations required for the method to converge, 
if it converges at all. The applicability of iterative methods 
is restricted since global convergence is not inherent to all 
initial matrices. A general condition for initialization is given 
by ||I — CC 0 II 2 < 1 or IIS 0 II 2 < 1- This condition ensures 
that the residual converges towards zero after each iteration. 

However, there are some conventional initialization methods 
which guarantee convergence. In lfl4l . theorem 2 shows that 
to find the inverse of a matrix C, the initialization Cq = 
aC^, where a satisfies 0 < a < — and denoted 

as the largest eigenvalue of the matrix A = C^C, ensures 
convergence. To reduce the complexity, following bound is 
used ifT^ 

^max — Supper — TTl f(A 1) 2 (13) 

where m = and ^ 

selected as a = 2/A„pper, which ensures convergence. In the 
next section, we propose a low complexity SD algorithm for 
large-antenna and massive MIMO systems that uses above ma¬ 
trix inversion methods to accurately estimate the transmitted 
signal vector. 

IV. Sphere Decoding using Iterative Matrix 
Inverse 

Now, let us investigate the advantages of using iterative 
matrix inverses for non-linear detectors, such as SD. Presently, 
there are two main versions of SD. The first is the Schnorr- 
Euchner enumeration Qa, m that updates the radius for 
SD adaptively, where after starting with an infinite radius, 
the search space shrinks with each good point until we get 
the optimal solution. In large/massive MIMO systems, such a 
technique would result in a huge decoding complexity. The 
other one is Fincke-Pohst algorithm based SD m, ini, 
which uses a fixed radius approach, and all the points that are 
inside the search space defined by the radius are compared for 
detecting the transmitted signal. This technique is extremely 
sensitive to the choice of the radius. It has been shown in 
the literature that both these approaches provide near ML 


Algorithm 1: Proposed SD Scheme 
Input ; y, H, H, k 

Output ; X 

Initialization i = K, cost = Vk, Ci = 0; 

[Q R] QR decomposition of H and z = Q^y; 

X ^ DFTS(z, R, H, cost, Ci, d, i); 

Function: DFTS(z, R, H, cost, Ci, i) 
for j ^ 1 to length{n) do 
I Cj = \zi — yxj € H; 

end 

Sort Cj’s in ascending order and keep only those symbols 
for which cj < {cost — c)); 
if Cj ^ {cost — Ci) then 
I return x, cost; 
else 

for u ^ 1 to length{c) do 

Xi — Xu, 

Cz Ci “t" Cu, 

if j = 1 then 

if costtemp < cost then 
cost Ci; 
return x, cost; 

end 

else 

Z - Z R,;^ZZ ^u. 

Extend the tree T for all H; 

[x, cost] DFTS(z, R, ft, cost,Ci, i — 1); 

end 

end 

end 


performance. In this section, we propose a mechanism to 
reduce the complexity of SD. 

Our SD algorithm combines both the strategies wherein we 
initialize with a BR computed using a low complexity iterative 
matrix inverse and also update the radius adaptively with every 
good point. The number of updates when using this algorithm 
would be significantly less, as the radius will be updated only 
when a new point is closer to the transmitted signal than ZF. 
Also, we are always guaranteed a solution as the ZF solution is 
always inside the searched domain. In Algorithm [T] we show 
the steps of the proposed SD scheme. 

A. Comparison of Bahai Radii Calculated through Approxi¬ 
mate and Exact Matrix Inverses 

Though iterative methods provide a good approximate in¬ 
verse, it is important to analyze the effect of approximation 
on the BR, as the choice of radius largely governs the 
complexity of SD. Interestingly, we show that the application 
of approximate matrix inversion methods also reduces the 
value of radius which leads to further savings in complexity. 

To prove this, let us define re as the BR computed through 
exact inverse and as the BR computed through the iterative 














Fig. 1. BR with Newton’s iterative method for a 16 x 16 system. 


method after k iterations. We know that re = limfe_>oorfe. 
Now, from the definition of BR ifTTl . ifT^ . we can write 

Te = ||R(Xzf -Xi)||, (14) 

rk = ||R(xzf - X2)||, (15) 

where Xi = C~^g, X 2 = C^g and R is obtained from the 
QR decomposition of H as H = QR. Let n denote the noise 
with respect to xzf, i-S-, 

y = HxzF + n. (16) 

Now, let us define the relation between transmitted vector x 
and detected ZF vector Xzf as 

X = Xzf + A, (17) 

where ||A|| denotes the magnitude of error in xzf- Since x 
and Xzf both belong to the same constellation, expectation of 
the difference between x and xzf would be zero. Substituting 
dnii in (fThl) . we get n = n — HA, and thus E(n) = 0. 

For a sufficient number of iterations, we can write the 
expected difference between the squares of the two radii in 
(fT^ and (fTsT i as 

nrl - rl] = 25R[E{n^HSfcC-ig}]. (18) 

We prove the above equation in Appendix I. 

We next show that the L.H.S. in (fTSl l decreases as the 
number of iterations increase. Using g = H^y in (fT^ . we 
can write 

E[n"HSfeC-ig] = E[n"HSfeC-m^(HxzF + n)] 

= E[n'f^HSfcXzF] + E[n^HSfeC”m^n]. (19) 

Also, for a given channel matrix H and received vector y, xzf 
would be a constant. Therefore, we can take vector xzf out 
of the first expectation term in ( fT9] l and it becomes 

E[n"HSfeXzF] = E[n^]HSfeXzF = 0, (20) 



Fig. 2. Bit error performance for the MMSE decoder in a massive MIMO 
system with N = 128, X = 8 for 16-QAM. 

and therefore, can be rewritten as 

E[n^HSfcC”ig] = E[n^HSfcC”m^n] 

= iVoTr(HSfeC”iH^), (21) 

where Tr(X) denotes the trace of matrix X. From ( fTSl l and 

(EB 

nrl-rl] = 27VoK[Tr(HSfeC-iH^)] 

= 2WoK[Tr(SfeC-^H^H)] 

= 2Wo5i[Tr(Sfe)]. (22) 

Similarly, for the radius obtained after fc +1 iterations, we get 

E[r2-r2+i] = 27VoK[Tr(Sfe+i)]. (23) 

It can be seen that for the residual matrix = I — C^C, 
we have Tr(Sfc) > 0. If the iterative methods used for matrix 
inversion converges to the exact inverse, it can be assumed that 
Tr(Sfe) > Tr(Sfe+i), as the elements of the residual matrix 
will tend towards zero as the number of iterations increase. 
Therefore, from equations (l22l i and (|2^ . it can be deduced 
that 

nrl-rUi] < nrl-rl] 

> E[rl] 

which means that, in general, BR after k iterations is smaller 
than the BR calculated after fc + 1 iterations. In Fig. [T] we 
use Newton’s iterative method for computing the approximate 
inverse and plot the BR for different iterations for a 16 x 16 
MIMO system. A monotonic rise in the value of BR with 
increasing iterations corroborates the above analysis. 

As Te = linifc_>oo rk, therefore 

E[rg] > E[r^] for finite k, 

i.e. the BR calculated using the exact inverse is larger than the 
BR calculated through an iterative method for all the iterations. 
Thus, as stated before, an approximate inverse can provide 
twofold savings in complexity. 

































































Fig. 3. Bit error performance for the ZF decoder in a massive MIMO system 
with TV = 128, X = 8 for 16-QAM. 



Fig. 4. Average number of computations for different SD schemes for a 
16 X 16 large MIMO system for 4-QAM. 

V. Simulation Results 

We first examine the performance of ZF and MMSE detec¬ 
tors for massive MIMO scenarios. Subsequently, we compare 
the simulation results for different SD methods that exist in 
the literature to the scheme proposed in Algortihm[T] In Fig.|2] 
and Fig. |2 we plot Bit-Error-Rates (BER) for MMSE and ZE 
decoders, respectively, for the cases when the matrix inverse 
is calculated exactly and using Newton’s iterative method. We 
calculate the approximate inverse for 3, 5 and 7 iterations. 
We see that for 3 and 5 iterations, the error performance 
in the case of MMSE is slightly away from the case when 
the exact inverse is used. However, increasing the number 
of iterations to 7 provides identical performance. Using more 
number of iterations would not result in any performance gain. 
Similarly, in the case of ZE decoding, performance improves 
with the number of iterations, and 7 iterations provides the 
same performance as the ZE decoder using the exact inverse. 

We also perform Monte Carlo simulations for BER and 
average number of computations for the three different SD 
schemes discussed above. The hrst two are adaptive radius 
(SE-SD) and hxed radius (FP-SD) algorithms respectively. 



Fig. 5. Bit error performance for different SD schemes for a 16 x 16 large 
MIMO system for 4-QAM. 



Fig. 6. Average number of computations for different SD schemes for a 
massive MIMO system with 32 base antennas and 8 users for 4-QAM. 

We compare these conventional schemes with the SD scheme 
proposed in Algorithm [T] in terms of performance and average 
number of computations required to hnd the solution. We use 
Newton’s iterative method with 7 iterations to calculate the 
approximate matrix inverse. In Fig.H] we compare the average 
number of computations required by the three schemes for 
a 16 X 16 system. It can be observed from the figure that 
the proposed SD scheme takes at least 35% less number of 
computations compared to the other two schemes. Also, from 
Fig. |5] we can deduce that there is no reduction in the quality 
of performance as all the three schemes give the same BER. 
In Eig. |6] and Eig. |7] we present similar numerical results for 
an V = 32 and K = 8 massive MIMO system. We again 
note that our SD scheme outperforms the conventional scheme 
while providing the same error performance. 

VI. Conclusion 

We have shown the advantages of using an approximate 
matrix inverse for detectors in large/massive MIMO systems. 
We obtained the maximum error which can be tolerated in the 
inverse to arrive at the same quantized ZE/MMSE solution. 





































































































Fig. 7. Bit error performance for different SD schemes for a massive MIMO 
system with 32 base antennas and 8 users for 4-QAM. 

Simulation results show that iterative inversion methods, used 
to calculate the ZF and MMSE solutions, reached the same 
performance as provided by the exact inverse for sufficient 
number of iterations. Extending the idea to complex detectors 
like SD, we show that the value of BR calculated using iter¬ 
ative methods is less than the BR obtained through the exact 
method. To this end, we proposed an adaptive SD scheme that 
uses BR as the initial radius. Simulation results show that the 
proposed SD scheme outperforms EP-SD and SE-SD in terms 
of complexity without any loss in performance. 

Appendix 

To prove (fT^ . we use the definition of Ve and rk from (fT4l i 
and (fTSl) so that 

rl-rl = ||R(xi-X zf)|P - ||R(x2-xzf)|P 

= ||RC-ig||2-||RCfcg||2 

+ 2{5R[(RxzF)^R(Cfe - C-i)g]} 

Now, using the fact that Ck = + Efc, we get 

rl-rl= 2K[(xzf - C-ig)"R^REfcg] - HRE^glp. (24) 

After using ( fTSl ) and taking expectations on both sides, we get 

E[rl - rl] = E[2K{(-C-iH"n)"R"REfcg} - ||REfcg|| = ]. (25) 

We will be neglecting the second term in R.H.S. of (|25]) citing 
the following assertion 

||REfeg||2 = (REfeg)^(REfeg) = g^EfR^REfeg. (26) 

Erom the orthogonal property of Q, we have R^R = 
H^H = C and therefore (l26l l becomes 

||REfcg||2 = g^EfCEfeg 
Using Efe = Cfc — C“^, we have 

IIREfcglp = g^(Cf-(C-i)")C(Cfe-C-i)g 
= g^(CfC^-I)(CfeC-I)C-ig 
= g^SfSfcXi, 


where = I — C^C is the residual matrix. Here, we have 
used the fact that C is a Hermitian matrix and (C“^)^ = 
C~^. Also, the first term in the R.H.S. of (l25T l can be written 
as 

E[23fi{(-C-iH^n)^R^REfcg}] = E[25ft(n^HEfcg)] 

= E[23?(n^HSfcC-ig)]. 

Eor sufficient number of iterations, would be very small and 
hence the term E[||REfcg^||] can be neglected when compared 
to first term in the R.H.S. of equation (l25T l. as the former is 
proportional to while the latter is proportional to Sk- 

Hence, (l25T l can be rewritten as 

nrl-rl] « E[23?(n^HSfeC-ig)] 
which proves (fTSl l for sufficient number of iterations k. 
References 

[1] E. Larsson, O. Edfors, F. Tufvesson, and T. Marzetta, “Massive MIMO 
for next generation wireless systems,” Communications Magazine, IEEE, 
vol. 52, no. 2, pp. 186-195, February 2014. 

[2] F. Boccardi, R. Heath, A. Lozano, T. Marzetta, and R Popovski, “Five 
disruptive technology directions for 5G,” Communications Magazine, 
IEEE, vol. 52, no. 2, pp. 74-80, February 2014. 

[3] L. Lu, G. Li, A. Swindlehurst, A. Ashikhmin, and R. Zhang, “An 
Overview of Massive MIMO: Benefits and Challenges,” Selected Topics 
in Signal Proc., IEEE Journal of, vol. 8, no. 5, pp. 742-758, Oct 2014. 

[4] J. Hoydis, S. ten Brink, and M. Debbah, “Massive MIMO in the UL/DL 
of Cellular Networks: How Many Antennas Do We Need?” Selected 
Areas in Communications, IEEE Journal on, vol. 31, no. 2, pp. 160- 
171, February 2013. 

[5] H. Li and V. Leung, “Low complexity zero-forcing beamforming for 
distributed massive MIMO systems in large public venues,” Comm, and 
Networks, Journal of, vol. 15, no. 4, pp. 370-382, Aug 2013. 

[6] A. K. Sah and A. K. Chaturvedi, “Reduced neighborhood search algo¬ 
rithms for low complexity detection in MIMO systems,” in GLOBECOM 
2015. IEEE International Symposium on, Dec 2015. 

[7] N. Srinidhi, T. Datta, A. Chockalingam, and B. Rajan, “Layered Tabu 
Seaiuh Algorithm for Large-MIMO Detection and a Lower Bound on 
ML Performance,” Communications, IEEE Transactions on, vol. 59, 
no. 11, pp. 2955-2963, November 2011. 

[8] X. Peng, W. Wu, J. Sun, and Y. Liu, “Sparsity-boosted detection for 
lai'ge MIMO systems,” Communications Letters, IEEE, vol. 19, no. 2, 
pp. 191-194, Feb 2015. 

[9] J. W. Choi and B. Shim, “New approach for massive MIMO detection 
using spai'se error recovery,” in Global Communications Conference 
(GLOBECOM), 2014 IEEE, Dec 2014, pp. 375^3759. 

[10] Y. Wang and H. Leib, “Sphere decoding for mimo systems with newton 
iterative matrix inversion,” Communications Letters, IEEE, vol. 17, no. 2, 
pp. 389-392, February 2013. 

[11] B. Hassibi and H. Vikalo, “On the sphere-decoding algorithm I. Ex¬ 
pected complexity,” Signal Processing, IEEE Transactions on, vol. 53, 
no. 8, pp. 2806-2818, Aug 2005. 

[12] E. Agrell, T. Eriksson, A. Vardy, and K. Zeger, “Closest point search in 
lattices,” Information Theory, IEEE Transactions on, vol. 48, no. 8, pp. 
2201-2214, Aug 2002. 

[13] F. Soleymani, “On a Fast Iterative method for Approximate Inverse of 
Matrices,” Communications of the Korean Mathematical Society, 2013. 

[14] V. Pan and R. Schreiber, “An improved newton iteration for the general¬ 
ized inverse of a matrix, with applications,” SIAM Journal on Scientific 
and Statistical Computing, vol. 12, no. 5, pp. 1109-1130, 1991. 

[15] T. Cui, S. Han, and C. Tellambura, “Probability-Distribution-Based Node 
Pruning for Sphere Decoding,” Vehicular Technology, IEEE Transactions 
on, vol. 62, no. 4, pp. 1586-1596, May 2013. 

[16] C. SchnoiT and M. Euchner, “Lattice basis reduction: Improved practical 
algorithms and solving subset sum problems,” Mathematical Program¬ 
ming, vol. 66, no. 1-3, pp. 181-199, 1994. 

[17] M. P. U. Fincke, “Improve Methods for Calculating Vectors of Short 
Length in a lattice, including a complexity analysis,” Math. Comp., 1985. 






















